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The  penetration  damage  created  in  carbon  composites  with  a range  of 
microstructures  has  been  characterized  for  both  quasi-static  and  dynamic 
conditions.  Several  quasi -static  properties  of  the  composites  have  also  been 
measured.  Analyses  of  the  quasi-static  penetration  have  been  developed,  which 
appear  to  present  a consistent  interpretation  of  the  damage  and  its  relation 
to  the  quasi-static  properties.  A preliminary  numerical  computation  of  the 
dynamic  damage  has  been  conducted,  with  encouraging  results.  Suggestions  for 
microstructural  modifications  that  should  ennance  the  damage  resistance  have 
been  made,  based  on  the  results  of  both  the  observations  and  analyses. 
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Studies  of  the  damage  in  carbon  composites  created  by  the  impact  of 
projectiles  have  identified  kinking  as  a primary  damage  mechanism.  The 

presence  of  kinks  below  and  adjacent  to  the  impact  crater  (Fig.  1)  apparently 
degrades  the  mechanical  strength  of  the  material,  and  permits  ready  ejection 
of  the  kinked  zone  upon  subsequent  impact.  The  characterization  and 
interpretation  of  the  microstructural  dependence  of  the  kink  zone  thus  has 
important  implications  for  the  design  of  composites  with  superior  impact 
damage  and  erosion  resistance. 

The  present  study  addresses  the  microstructural  issue  by  examining 
the  projectile  damage  created  in  carbon  composites  with  a range  of 
microstructures.  The  complexity  of  impact  damage  problems  is  recognized  and 
has  been  described  in  some  detail  for  carbon-cargon  composites. ^ In 
consequence,  following  an  approach  adopted  for  the  analysis  of  Impact  damage 
in  Isotropic  brittle  materials,  a series  of  carefully  selected  quasi-static 
studies  are  firstly  performed  and  analyzed,  to  provide  a qualitative 
perspective  from  which  the  impact  problem  can  be  viewed.  The  quasi-static 
studies  Include  compression  tests  designed  to  examine  kink  initiation, 
compliance  tests  to  assess  the  constituent  properties,  and  penetration  tests 
to  determine  the  dynamic  penetration  resistance. 

The  quasi-static  studies  are  followed  by  projectile  impact  tests, 
which  characterize  the  extent  of  kinking  and  Its  dependence  on 
microstructure.  A sample  numerical  calculation  is  then  presented  which 
Indicates  an  approach  for  analyzing  the  Impact  damage  phenomenon  in  detail. 
Finally,  the  implications  of  the  study  for  both  the  development  of  simple 
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2.0  MATERIALS 

Carbon  composites  of  the  type  considered  in  this  paper  consist  of 
bundles  of  carbon  fibers,  arranged  in  an  orthogonal  Interpenetrating  network, 
disposed  in  a matrix  of  graphite.  The  graphite  thus  exists  as  interfiber 
material  within  the  bundles  and  as  Isolated  interstices  between  bundles. 

The  graphite  matrix  invariably  contains  porosity  and  often  the  interstices  are 
separated  from  one  (or  more)  of  the  adjacent  bundles.^  A proportion  of 
the  individual  fibers  within  the  bundles  are  fractured,  at  statistically 
distributed  locations;  frequently,  regions  of  considerable  fiber  distortion 
are  in  evidence.^  The  role  of  these  microstructural  inhomogeneities  is  an 
Important  theme  of  this  paper. 

The  characteristics  of  the  materials  selected  for  detailed  study  are 
summarized  in  Table  I.  The  rationale  for  selecting  these  materials  was  based 
on  the  premise  that  damage  trends  could  be  most  effectively  deduced  if  the 
microstructural  (bundle  and  interstice)  dimensions  and  the  properties  of  the 
constituents  (fibers  and  matrix)  could  be  independently  varied.  Consequently, 
materials  A,  B,  and  C are  observed  to  have  similar  relative  bundle  areas,  but 
different  constituents;  while  materials  C,  0,  and  E have  similar  constituents 
but  different  bundle  areas. 
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3.0  QUASI -STATIC  TESTS 

Previous  studies^*3 4)  have  indicated  that  kinking  can  occur  when 
the  fiber  bundles  are  subjected  to  longitudinal  compression,  lateral  shear  or 
lateral  compression:  each  mode  being  prevalent  in  impacted  samples.  The 
available  analyses  of  these  kink  formation  processes, although  capable 
of  rationalizing  certain  aspects  of  kink  formation,  are  still  too  rudimentary 
to  predict  kink  initiation  conditons.  An  experimental  measurement  is  thus 
required.  A constrained  compression  test  designed  to  provide  this  information 
is  described  in  this  section.  The  interpretation  of  the  test  results  invokes 
analyses  of  compressive  bundle  deformation,  described  in  Appendix  I,  and  of 
kink  initiation,  described  in  Appendix  II. 

The  penetration  resistance  afforded  by  a particular  microstructure  is 
related  in  some  complex  fashion  to  the  kinking  process,  and  other  mechanisms 
of  deformation.  A measure  of  the  penetration  resistance  Is  obtained  by 
performing  penetration  tests,  and  some  of  the  controlling  parameters  are  eluci- 
dated by  examining  the  damage  surrounding  the  penetration  zone.  A damage 
analysis  is  also  undertaken,  to  provide  insights  pertinent  to  the  Impact 
damage  problem.  The  analysis  utilizes  approximate  relations  for  the 
deformation  that  accompanies  penetration,  developed  in  Appendix  III. 

3.1  CONSTITUENT  MODULI 

The  elastic  and  plastic  properties  of  the  mlcrostructural 
constituents  are  important  parameters  for  all  Impact  damage  models.  A series 
of  microcompression  tests  have  been  designed  to  obtain  approximate  values  of 
these  Important  properties.  The  values  are  necessarily  approximate  in 
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character  because  each  property  is  derived  from  uniaxial  force/displacement 
curves  (compliance  curves)  by  assuming  elastic  isotropy.  More  precise 
measurements  of  the  individual  constituent  properties  will  ultimately  be 
needed.  However,  this  is  a nontrivial  problem;  which  is  not  attempted  in  the 
present  study,  since  the  current  objective  is  to  obtain  a measure  of  the 
relative  trends  in  these  properties  with  microstructure. 

An  in-sltu  compliance  measurement  technique  has  been  described 
previously, ^ and  provides  values  for  the  constituent  properties  summarized 
in  Table  II.  In  this  study,  the  in-situ  tests  were  augmented  by  individual 
compressive  compliance  tests  on  small  sections  of  each  mlcrostructural 
constituent.  The  constituents  were  Isolated  using  a precision  dicing  saw  and 
compression  tests  conducted  on  a small  Instron  machine  (with  0.1  kg  load 
cell),  using  an  LVTD  clip  gauge  device  to  monitor  the  displacement  of  the 
sample.  The  results  are  sunmarlzed  in  Table  II.  It  is  apparent  that  the  two 
methods  yield  comparable  results. 


3.2  CONSTRAINED  COMPRESSION  TESTS 

A quasi-static  technique  for  measuring  the  onset  conditions  for 
kinking  under  longitudinal  compression  has  been  devised.  The  technique 
(Fig.  2a)  involves  the  application  of  equal  steady  state  compressive 
contralnts  in  two  orthogonal  directions,  with  a monotonlcally  Increasing 
compression  applied  in  the  third  orthogonal  direction.  The  use  of  constraints 
is  Important  because  this  ensures  that  the  damage,  when  it  Initiates,  occurs 
as  kinks  (Fig.  3):  rather  than  bundle  fracture,  fiber /matrix  debonding,  etc., 
which  frequently  occur  in  the  absence  of  constraint.  The  apparatus  used  for 
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this  experiment  is  depicted  in  Fig.  2b.  The  constraints  are  applied  through 
mechanical  screw  devices,  and  the  magnitudes  of  the  forces  applied  to  the 
composite  are  measured  using  strain  gauges  attached  to  the  compression 
blocks.  The  orthogonal  compression  is  applied,  through  a self  aligning 
cylindrical  plunger,  by  a rigid  constant  displacement  rate  device  (e.g.,  an 
Instron  testing  machine).  The  load  and  displacement  are  recorded  using 
conventional  load  cell  and  displacement  gauge  (LVDT)  instrumentation. 

A typical  load  displacement  behavior  obtained  under  constrained 
compression  is  shown  in  Fig.  4.  Initially,  the  load  increases  linearly  with 
displacement.  Subsequently,  there  is  a small  load  decrement  (yield  point), 
followed  by  a nonlinear  Increase  in  load  (to  a load  maximum).  Sectioning  of 
samples  tested  to  various  locations  on  the  load/deflection  curve  indicates 
that  there  is  no  detectable  damage  in  the  linear  region.  However,  immediately 
above  the  load  discontinuity,  a kink  Is  invariably  detected  in  one  of  the 
longitudinal  bundles  (Fig.  3).  It  is  presumed,  therefore,  that  this 
discontinuity  Is  associated  with  the  onset  of  longitudinal  bundle  kinking.  At 
larger  loads,  well  into  the  nonlinear  region  (Fig.  4),  complex  kinking  of  the 
lateral  bundles  is  first  discerned. 

Critical  values  of  the  kink  Initiation  conditions  are  obtained  from 
the  load  discontinuity,  for  each  of  the  materials  listed  In  Table  I.  The 
critical  kink  initiation  stresses  are  obtained  from  the  applied  strain  at 
yield  and  the  approximate  moduli  of  the  bundles  and  matrix  (Table  II).  The 
critical  stresses  are  plotted  on  Fig.  5.  For  the  three  materials  (A,B,C)  with 
relatively  isotropic  bundle  dimensions,  the  critical  stresses <r*J  are 

w 

independent  of  the  constraint.  However,  Increases  with  Increase  in 
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constraint  for  the  two  materials  with  a predominance  of  longitudinal  bundles 
(D,E).  Sectioning  of  these  samples  after  testing  indicated  that  the  lateral 
bundles  in  the  latter  two  samples  had  been  extensively  kinked  by  the 
constraining  loads  (thereby,  enabling  the  constraints  to  be  transmitted  to  the 
longitudinal  bundles). 

The  results  obtained  at  the  minimum  constraint  (~15  MPa)  are 
replotted  in  Fig.  6 in  terms  of  the  cumulative  kink  initiation  probability, 
#(cr'),  obtained  from  order  statistics. 

Kink  formation  under  axial  compression  is  associated  with  the 
presence  of  inhomogeneities  within  the  bundles  (Appendix  I).  The  involvement 
of  inhomogeneities  has  statistical  implications  for  the  kink  initiation 
stress,  reminiscent  of  the  statistics  of  brittle  fracture. Specifically, 
for  a spatially  random  distribution  of  inhomogeneities,  the  probability,  p,  of 
kink  initiation  in  a small  element  of  bundle,  length  Ai  and  cross  section 
(s 

f 9(Sk)dSk 
0 

where  g ( Sk ) dS^  is  the  number  of  inhomogeneities  in  unit  volume  of  bundle 
that  produce  kinks  at  applied  axial  compressive  stresses  between  and 
Sk  + dS^,  and  cr  is  the  level  of  the  applied  compression.  Then,  if  there 
Is  no  Interaction  between  Inhomogeneities,  the  probability  of  kink  Initiation 
Is  bundles  of  total  length  1 can  be  obtained  from  the  product  of  the  survival 
probabilities  of  each  element,  as: 
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#(o£) 


- p)  ■*  1 - exp 


r ar°r 

dl  /9<Sk>dSk 

Ji  0 J 


(2) 


For  the  present  tests,  the  applied  stress  Is  uniform  and  Eq.  (2)  reduces  to: 


In 


1 -#(«£) 


V J 9(S|c)dS|{ 


(3) 


The  function  g($K)  can  also  be  expressed  in  terms  of  £,  obtained  by  direct 
inversion  of  Eq.  (3),  as: 

g(Sk)  - (d^/dSk)/Abj€  (4) 


The  values  of  g($k),  obtained  *rom  a polynomial  fit  of  the  data  in  Fig.  6, 
are  plotted  in  Fig.  7.  It  is  apparent  from  Fig.  7 that  a suitable  function 
for  g(Sk)  is; 


9(Sk)  - \(sk/s*)k 


(5) 


where  k is  a shape  parameter,  S*  is  a scale  parameter,  and  \ is  a constant. 

The  specific  values  of  k,  S*  and  \ pertinent  to  each  material  are  indicated  on 
Fig.  7,  These  values  of  g(Sk)  can  hereafter  be  inserted  into  Eq.  (2)  to 
obtain  kink  initiation  probabilities  for  inhomogeneous  axial  deformation. 


3.3  PENETRATION  TESTS 

An  Instrumented  drop  tower  was  used  for  the  penetration  tests^ 
(this  test  might  be  considered  dynamic  because  it  involves  penetration 

7 

C/Q726A/ebs 


Rockwell  International 

Science  Center 


SC5076.5FR 

velocities  of  a few  meters  per  second,  four  orders  of  magnitude  greater  than 
conventional  static  loading  rates).  By  directly  Instrumenting  the  penetrator, 
reasonably  accurate  simultaneous  measurements  of  the  load  and  deflection  could 
be  obtained.  A deflection  gauge  system  was  used  for  penetration 
measurements.  Two  Biomatlon  Model  802  Transient  Records  were  used  to  record 
load  and  deflection  signals.  The  records  were  synchronized  by  a 
Tetronix  FG  501  Function  Generator.  The  Information  was  digitized,  fed  Into  a 
Hewlett-Packard  Model  9830A  computer,  and  stored  on  magnetic  tape.  The  raw 
data  could  then  be  plotted  or  further  processed.  The  final  result  is  a plot 
of  penetration  pressure,  p,  versus  depth  of  penetration  q/rj  (normalized  to 
the  indenter  radius  r^).  A curve  fitting  procedure  defined  by 

P - P0  (q/2r j ) n (6) 

was  then  applied  to  the  data.  The  coefficient  pQ  is  the  pressure  at  one 
diameter  of  penetration,  while  the  exponent  n is  the  slope  of  the 
pressure/penotration  curve  (related  to  the  work  hardening  exponent).  The 
parameters  pQ  and  n are  determined  by  a least  mean  squares  fit.  Penetration 
pressures,  pQ,  and  work  hardening  exponents  are  obtained  for  each  material 
are.  summarized  in  Table  III. 

The  damage  In  the  penetration  samples  was  characterized  optically, 
after  sectioning  through  the  penetration  zone  and  polishing.  The  damage 
obtained  with  a relatively  large  (3  mm  diameter)  penetrator  Is  summarized  In 
Figs.  8 and  9.  The  damage  In  the  longtudinal  bundles  has  several  important 
features  (Fig.  8).  The  kinking  In  the  bundles  beneath  the  impact  crater 
(bundles  C and  D,  and  B,  E,  F below  lateral  bundle  3)  appears  to  exhibit 
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little  spatial  correlation,  either  along  the  bundles  or  between  adjacent 
bundles.  Kink  pairs  of  opposite  sign  appear  to  be  the  primary  mechanism  for 
retaining  the  Initial  bundle  alignment.  The  kinks  at  the  extremity  of  the 
damage  zone  generally  form  at  inhomogeneities,  as  Indicated  on  the  figure. 

The  kinking  in  the  bundles  beside  the  crater  (bundles  B and  F)  tend  to  exhibit 
a uniform  spatial  separation,  with  the  kinks  usually  initiating  near  the 
corners  of  an  out-of-plane  lateral  bundle  (e.g.,  lb,  2b,  If,  2f,  and  3f).  The 
kink  separation  is  thus  approximately  equal  to  the  diameter  of  a lateral 
bundle  plus  a matrix  Interstice.  Also,  the  kinks  tend  to  conform  to  a single 
orientation.  Finally,  it  is  noted  that  kinks  occasionally  occur  in  bundles 
that  are  not  directly  subject  to  Interaction  with  the  penetrator  (e.g., 
bundle  A). 

The  damage  to  the  in-plane  lateral  bundles  also  has  important 
implications  (Fig.  9).  The  lateral  extent  of  the  kinks  in  the  bundles  beside 
the  crater  is  restricted  to  less  than  a few  bundle  diameters.  Usually,  these 
kinks  are  highly  correlated,  initiating  from  the  corners  of  out-of -plane 
lateral  fiber  bundles  along  a single  plane.  It  is  instructive  to  regard  these 
kinks  as  being  induced  by  a shear  displacement  as  depicted  by  the  vectors  in 
Fig.  9.  Complex  kinks  usually  appear  directly  beneath  the  impact  center  and 
occasional  fractures  occur  In  the  same  vicinity  (Fig.  9). 

The  depth  of  the  kink  zone  in  the  longitudinal  bundles  depends  on 
both  the  diameter  of  the  penetrator  end  the  penetration  rate.  The  depths 
below  the  crater  measured  for  a 1 mm  diameter  penetrator  are  summarized  In 
Table  III  and  on  Fig.  16. 
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3.4.1  Compression  Tests 

The  kink  Initiation  stresses  obtained  in  the  compression  tests  can  be 
rationalized  in  terms  of  the  enhanced  curvature  produced  by  an  applied 
compression  at  local  Inhomogeneities  (Appendix  I),  coupled  with  the  critical 
kink  formation  conditions  (Appendix  II).  The  fiber  curvatures  that  result 
from  an  applied  compression  depend  most  prominently  on  the  initial  dimensions 
(y,  1 ) of  the  inhomogeneity  (Appendix  I),  although  the  matrix  yield  strength 
and  lateral  bundle  modulus  can  also  be  significant.  The  fiber  curvatures 
generate  local  tensile  stresses,  which  are  the  driving  force  for  kink 
initiation.  The  magnitude  of  the  tensile  stress  required  to  initiate  kinking 
is  typically  influenced  by  the  matrix  yield  strength  (Appendix  II)  with  the 
fiber  fracture  strength  and  modulus  also  being  important  in  certain  regimes. 
The  critical  applied  compression  for  kink  formation  should  thus  be  primarily 
Influenced  by  the  inhomogeneity  dimensions  and  the  matrix  yield  strength,  with 
the  constituent  moduli  and  fiber  fracture  strength  as  possible  alternatives. 
The  trends  in  the  kink  initiation  stress  are  broadly  Interpreted  in  terms  of 
these  influences. 

The  relatively  Invariant  kink  formation  stress  <rj$  (at  low 
constraints)  obtained  for  the  three  materials  containing  the  same  fiber  and 
matrix  constituents  (C,D,E)  Indicates  that  there  are  no  unexpected  influences 
of  the  bundle  dimensions.  The  elevation  ofcrjj  by  lateral  constraint 
(Fig.  5),  when  this  constraint  Is  transmitted  to  the  longitudinal  bundles,  can 
be  Interpreted  as  resulting  from  a suppression  of  matrix  yield  (i.e.,  through 
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a Von  Mlses  criterion,  Eq.  (12))  or  perhaps  a decrease  In  the  Initial 
amplitude  y of  inhomogene  It les.  This  would  suggest  a criterion  for  kink 
formation  that  depends  on  the  difference  bo  between  the  maximum  a-  and 
minimum  <r^n  compressive  stress  (i.e.,  the  shear  stress).  Additional  work 
Is  needed  to  establish  the  detailed  nature  of  this  criterion;  although  studies 
on  uniaxial  carbon/epoxy  composites^  infer  a direct  proprotionallty 
UcraJ$).  For  purposes  of  further  analysis  it  is  thus  assumed  that 
kinking  initiates  when  <r£  is  equal  to  the  measured  kink  initiation  stress 
at  zero  constraint  (Figs.  5,  6). 

The  strong  effects  of  the  fiber  and  matrix  properties  on  the  kink 
initiation  stress  (materials  A,  B,  C)  are  of  great  significance.  Inspection 
of  the  material  microstructures  (cf.  Figs.  11  and  12)  indicates  that  the 
fibers  in  materials  A and  B exhibit  more  extensive  perturbations  from 
linearity  than  in  material  C,  since  only  in  the  latter  can  individual  fibers 
be  traced  for  significant  distances  on  a planar  section.  Also,  the  latter  has 
a much  reduced  density  of  voids,  both  within  the  bundles  and  within  the 
interstices,  accounting  for  its  greater  yield  strength  (Table  II).  The 
superior  material  (C)  is  thus  characterized  by  a reduced  size  distribution  of 
inhomogeneities  (smaller  y)  and  a larger  matrix  yield  strength.  Qualitative 
conformance  with  the  expectations  of  the  simple  kink  formation  analysis 
described  in  Appendices  I and  II  is  thus  apparent  (a  more  quantitative 
interpretation  is  not  warranted  at  this  stage).  These  results  indicate  that 
further  improvements  in  kink  resistance  should  accompany  a reduction  in  the 
size  of  the  larger  intrabundle  inhomogeneities  (i.e.,  void  regions  and  regions 
of  initial  curvature)  and  an  enhancement  of  the  matrix  yield  strength  (e.g.. 
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reduction  In  porosity,  substitution  or  addition  of  high  strength  matrix 
constituents) . 

3.4.2  Penetration  Tests 

The  penetration  damage  Involves  additional  considerations.  The  kinks 
in  the  axial  longitudinal  bundles  have  both  in-plane  and  out-of -plane 
characteristics:  the  in-plane  kinks  usually  occur  as  pairs  of  opposite  sign, 
allowing  the  lateral  displacement  of  the  bundle  to  remain  Invariant.  The  kink 
locations  within  the  bundle  do  not  relate  to  the  out-of -plane  orthogonal 
bundles,  indicating  that  there  is  little  damage  enhancement  associated  with 
these  features.  The  zor.s  directly  beneath  the  crater  Is  intensely  kinked;  the 
kinks  within  the  zone  do  not  appear  to  initiate  at  obvious  inhomogeneities. 
However,  toward  the  limit  of  penetration  of  the  kink  damage,  the  kinks  are 
almost  invariably  associated  with  pre-existent  bundle  inhomogeneities 
(especially  voids  and/or  bundle  curvature).  These  kink  characteristics  are 
consistent  with  the  damage  anticipated  for  axial  compressive  deformation  of  a 
bundle,  as  noted  above.  The  remote  kinks  occur  in  regions  of  relatively  low 
deformation  (Appendix  III)  and  must  seek  a severe  inhomogeneity  to  satisfy  the 
kink  formation  requirements;  while  the  kinks  close  to  the  projectile  Interface 
develop  in  regions  of  much  larger  stress,  and  can  thus  Initiate  at  more 
abundant  small  Inhomogeneities. 

The  interrelation  between  the  extent  of  kinking,  z,  the  penetration 
pressure  and  the  kink  formation  probability  can  be  found  directly  from 
Eqs.  (2),  (5),  and  (A30).  If  it  Is  assumed  that  the  remote  kinks  are 
Influenced  by  the  remote  elastic  stress  field  (Eq.  A30b),  the  probability* 
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that  kink  initiation  will  occur  In  an  axial  bundle,  outside  a zone  of  depth  z, 
is 


[i  -♦]  ■ \ 


'AT 

z 


n 


/p  (a/z)2 


(S/S*)KdSdz 


(7) 


Integration  and  rearrangement  then  gives; 


XAb(/3*P0a2 

jk+1 

(S*)k  (k+1)  (2k+l)  ( 

- In  (l4)) 

l/(2k+l) 


(8) 


By  choosing  a suitably  small  value  for  #(e.g.,  0.01),  z becomes  the  extent  of 
the  damage  zone.  It  is  Instructive  to  obtain  ratios  of  z at  constant  * for 
the  three  materials  with  comparable  bundle  dimension  (A,B,C),  by  inserting  the 
measured  kink  formation  parameters  (Fig.  7)  and  the  penetration  pressures 
(Table  III)  into  Eq.  (8).  The  resultant  ratios,  referred  to  material  C,  are 
indicated  in  Table  III.  The  correlation  with  the  measured  values  (obtained  at 
low  penetration  rates)  Is  quite  good,  indicating  the  self-consistency  of  the 
measurement  techniques. 

The  damage  in  the  lateral  bundles  derives  from  both  the  non-axial 
deflection  and  the  axial  deformation  (e.g.,  Eq.  A42).  The  deflection  stresses 
increase  as  the  projectile  Interface  is  approached  (x  decreases).  Yet,  the 
only  kinking  that  occurs  beneath  the  penetration  zone  (x  - 0)  is  the  complex 
kinking  associated  with  the  compression  between  orthogonal  lateral  bundles 
(Fig.  9 ).  Also,  fractures  are  occasionally  observed  in  this  zone  (Fig.  9). 
The  quasi-static  deformation  field  (Fig.  23)  indicates  axial  tension  in  this 
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region;  Implying  that  kink  formation  is  suppressed  by  axial  tension. 
Presumably,  the  tension  tends  to  reduce  the  severity  of  pre-existent 
Inhomogeneities  (reduces  y in  Eq.  A19)  and  prevents  kink  instability 
conditions  from  being  attained,  even  though  the  non-axial  bundle  deflections 
are  relatively  large.  In  principle,  the  incidence  of  fracture  In  the  various 
materials  could  be  ascertained  by  equating  the  bundle  stresses  In  the  tensile 
zone  (e.g.,  Eq.  A18)  to  the  fiber  fracture  probability.  Fracture  data  are  not 
presently  available  for  each  material,  and  the  comparison  cannot,  therefore, 
be  conducted  at  this  juncture. 

The  kinking  that  occurs  in  the  zone  of  axial  compression  (x/z  ^ 1A/2) 
tends  to  concentrate  at  orthogonal  lateral  oundles  (Fig.  9).  The  local 
deflections  at  these  bundles  must,  therefore,  exert  an  important  influence  on 
kink  initiation.  This  can  be  qualitatively  rationalized  by  noting  that  the 
curvature  due  to  the  local  deflection  (Eq.  A40)  decreases  less  rapidly  with 
distance  x from  the  penetration  axis  than  the  curvature  due  to  the  general 
deflection  (Eq.  A35);  Indicating  that  the  former  should  dominate  at 
appreciable  values  of  x.  The  role  of  the  axial  compression  Is  to  preserve, 
and  perhaps  accentuate,  the  Influence  of  local  inhomogeneities  on  the 
formation  of  a kink  Instability.  Assuming  that  axial  compression  only 
preserves  the  inhomogeneities,  inserting  the  measured  kink  initiation  stresses 
(Fig.  5),  penetration  pressures  (Table  III)  and  projectile  radii  (500  urn)  Into 
Eq.  A42,  the  lateral  extent  x of  kinking  In  each  material  can  be  estimated 
(E^/<  E > Is  taken  to  be  102,  Table  II).  The  analysis  indeed  confirms 
that  the  local  deflection  term  Is  dominant,  and  predicts  values  of  x * 2 mm, 
reasonably  consistent  with  the  present  observations  (Fig.  9). 
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The  kinking  In  the  non-axial  longitudinal  bundles  Is  expected  to 
differ  from  that  in  the  lateral  bundles  because  there  is  no  axial  tension  near 
the  surface  to  restrict  the  formation  of  kinks  at  locations  of  maximum 
curvature  (z  ~ 0).  Intense  kinking  can  thus  be  anticipated  near  the  surface, 
In  accord  with  the  observations  (Fig.  8). 

4.0  IMPACT  TESTS 


4.1  TEST  TECHNIQUES 

The  method  of  particle  impact  testing  is  an  exploding  foil 
apparatus^  that  permits  small  spherical  projectiles  to  be  propelled  toward 
a target  at  very  high  velocities  (up  to  ~ 7,000  ms-1)  with  excellent 
accuracy.  This  apparatus  was  used  to  propel  1 mm  diameter  glass  spheres  onto 
carbon  composite  targets  of  the  type  listed  in  Table  I.  The  targets  were 
typically  cylindrical  in  shape,  2.5  cm  in  diameter  and  ~3  cm  thick. 

After  testing,  damage  overviews  were  obtained  in  the  scanning 
electron  microscope.  Then,  the  targets  were  vacuum  impregnated  with  polymer, 
and  sections  through  the  crater  suitable  for  subsequent  optical  microscopy, 
were  prepared  by  polishing  with  slurries  of  fine  alumina  powders.  The 
subsurface  damage  on  successive  sections  was  then  pictorial ly  recorded,  using 
partially  polarized  light  microscopy. 

4.2  DAMAGE  CHARACTERIZATION 

The  damage  in  all  cases  consisted  of  a zone  of  material  removal  (the 
crater),  bundle  kinks  and  occasionally,  bundle  fractures.  Typical  craters 
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(Fig.  10)  consist^)  of  a central  deep  penetration  zone  and  a shallower 
outer  ledge.  The  dimensions  of  the  crater  vary  systematically  with  the  Impact 
velocity  (Fig.  10).  There  is  also  a significant  dependence  on  the  bundle 
dimensions^  but  no  discernable  effect  of  the  fiber/matrix  properties. 

The  kink  damage  surrounding  the  crater  exhibits  the  general  features 
depicted  in  Figs.  1,11-15.  The  effects  of  Impact  velocity  and  microstructure 
on  the  maximum  depth,  z,  of  kinking  in  the  longitudinal  bundles  (Figs.  11-14) 
are  summarized  in  Fig.  16.  The  kink  depth  increases  monotonically  with 
velocity,  depends  prominently  on  the  fiber/matrix  properties  and  marginally  on 
the  relative  bundle  dimensions.  An  important  observation,  comparable  to  that 
for  the  penetration  test,  is  that  the  longitudinal  bundle  kinks  appear  to  form 
homogeneously  immediately  beneath  the  crater  (i.s.,  their  density  is  high  and 
they  completely  traverse  the  bundle);  but  are  almost  invariably  associated 
with  microscopic  inhomogeneities  at  their  maximum  depth  of  penetration 
(Figs.  14a  and  b).  The  inhomogeneities  that  are  most  frequently  associated 
with  these  kinks  are:  missing  fiber  segments,  severly  distorted  fibers,  major 
fiber/matrix  disbonds,  matrix  voids.  An  unusual  mode  of  kinking  (Fig.  14a)  is 
also  observed  occasionally  in  materials  A and  B,  occurring  at  the  corners  of 
out-of-plane  lateral  bundles. 

The  corresponding  trends  in  the  maximum  lateral  extension,  x,  of 
kinks  in  the  lateral  bundles  (Fig.  15)  are  summarized  in  Fig.  17.  The  extent 
of  this  mode  of  kinking  is  very  sensitive  to  the  impact  velocity  and 
microstructure;  substantial  effects  of  velocity,  bundle  dimensions  and 
fiber /matrix  properties  are  apparent.  Complex  kinking  and  fracture  of  the 
lateral  bundles  is  also  apparent  beneath  the  penetration  axis  (Fig.  15). 
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4.3  DAMAGE  INTERPRETATION 

The  strong  effects  of  velocity  on  the  relative  extent  of  kinking  in 
both  the  axial  longitudinal  and  lateral  bundles  indicates  that  these  key  modes 
of  damage  cannot  be  interpreted  by  a simple  extension  of  the  quasi -static 
analysis  described  previously  (Section  3.4).  A fully  dynamic  analysis  is 
required.  The  complexity  of  the  dynamic  problem  requires  that  numerical 
techniques  be  developed  for  detailed  damage  analysis  and  interpretation.  An 
approach  is  described  In  the  following  section  (Section  5). 

A qualitative  interpretation  of  the  damage  in  terms  of  available 
(isotropic)  dynamic  solutions  can,  however,  provide  some  useful  Insights.  The 
ratios  of  the  damage  depths  z in  the  axial  longitudinal  bundles  , for 
materials  A and  B versus  material  C,  are  greater  at  high  Impact  velocities 
than  the  equivalent  quasi-static  values  (Fig.  16).  The  most  self-consistent 
interpretation  of  this  observation  is  that  the  spatial  dependence  of  the 
dynamic  stress  is  less  than  the  quasi -static  dependence,  providing  a greater 
dependence  of  the  damage  depth  on  the  kink  formation  stress  (see 
Section  3.4).  Such  behavior  Is  anticipated  by  the  solutions  for  a pressurized 
cavity. (Also,  the  high  modulus  fibers  could  lead  to  a partial 
wave-guide  phenomenon,  further  reducing  the  geometric  attenuation  typically 
attributed  to  spherical  waves.)  Specifically,  if  it  is  assumed  that  the  ratio 
of  dynamic  penetration  pressures  Is  similar  throughout  to  the  values  obtained 
in  the  penetration  tests  (Table  III),  a generalized  form  of  Eq.  (8) 
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where  n is  the  spatial  exponent  of  the  stress  field  (oj,a(a/z)n),  Indicates 
that  the  observed  z ratios  at  high  velocities  ($-4,000  ms’*)  are  consistent 
with  a linearly  decreasing  stress  field  (i.e.,  cr  « 1/z).  The  alternative, 
but  less  plausible,  interpretation  of  the  dynamic  influence  would  require  that 
the  relative  dynamic  penetration  pressures  be  substantially  different  from 
those  obtained  quasi-statically. 

The  even  stronger  velocity  dependence  of  kinking  in  the  lateral 
bundles  (Fig.  17)  suggests  that  additional  dynamic  phenomena  are  involved. 
Available  dynamic  solutions  indicate  that  the  dynamic  stress  field  for  a rigid 
projectile  and  a plastically  deforming  target  are  relatively  uniform  about  the 
center  of  impact. The  tendency  for  kinking  in  the  lateral  bundles  to  be 
comparable  to  that  in  the  axial  longitudinal  bundles  might,  therefore,  be 
anticipated.  The  observed  greater  extent  of  the  lateral  bundle  kinking  is 
thus,  presumably,  associated  with  the  non-axlal  displacement  of  the  lateral 
bundles,  coupled  with  the  dynamic  axial  compression.  The  observed  in-plane 
morphology  of  the  kinks  in  the  lateral  bundles^  (Fig.  15)  is  consistent 
with  this  Interpretation. 

Finally,  the  consistent  trends  in  the  effect  of  microstructure  on  the 
three  measured  properties--the  extent  of  the  kink  zone,  the  yield  strength  in 
constrained  compression,  and  the  penetration  pressure— is  re-emphasized. 
Further  study  Is,  of  course,  needed  to  establish  the  quantitative 
Interdependence  of  these  properties  (I.e.,  a dynamic  equivalent  of  Eq.  8); 
but,  meanwhile,  these  relatively  simple  quasi-static  tests  should  provide  a 
measure  of  a material's  kink  resistance,  suitable  for  Inclusion  in  materials 
development  studies. 
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5.0  NUMERICAL  CALCULATION  OF  KINKING 


A preliminary  attempt  at  the  numerical  calculation  of  kinking  In  the 
longitudinal  bundles  has  been  undertaken  using  an  Eulerian  finite  difference 
method  (the  HELP  code).^2^  The  calculation  was  conducted  using  the 
axisymmetric  configuration  shown  in  Fig.  18;  in  which  a 500  Mm  radius  glass 
sphere  was  directed  to  impact  the  axis  of  a cylindrical  rod,  having  the 
average  properties  of  a fiber  bundle  (Table  II).  The  surrounding  material  was 
taken  to  be  isotropic,  having  the  average  properties  of  both  the  graphite 
matrix  and  the  fiber  bundles  (see  Table  II). 

5.1  COMPUTATIONAL  SCHEME 

The  computational  grid  was  chosen  to  have  the  finest  resolution 
(Fig.  18b)  adjacent  to  the  bundle  axis,  and  at  a depth  where  experiment 
(Fig.  16)  had  Indicated  that  kink  damage  could  be  anticipated.  The  equation 
of  state  used  to  describe  the  material  response  was: 

p - Af  + B^2  (10) 

where  p Is  the  hydrostatic  pressure,  A and  B are  the  constants  listed  In 
Table  IV:  ( is  the  volumetric  strain,  which  is  related  to  the  density  change 
by; 

f-P/P0~l  (11) 

where  PQ  Is  the  Initial  density.  The  anisotropic  constitutive  model  applied 
to  the  bundle  has  been  documented  elsewhere. ^12^  The  glass  and  matrix 
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materials  were  modeled,  respectively,  as  isotropic,  linear  strain  hardening 
and  elastic/perfectly  plastic.  In  both  cases,  a von  Mises  yield  criterion  was 
assumed, 


Vu 


x crw 


(12) 


where  S^.  is  the  stress  deviator  tensor  and  <r^  is  the  plastic  flow  stress 
in  shear  given  by; 

°y  * l°y°  + <ryf  ) (1  * U/Um}  (13) 

<r°  is  the  yield  stress  at  zero  pressure,  a-'  is  a measure  of  the 
* </ 

pressure  dependence  of  the  yield  stress;  U,  the  specific  internal  energy,  and 
Um,  the  energy  at  melting,  are  introduced  to  provide  an  approximate  measure 
of  the  temperature  dependence  of  the  yield  stress.  The  values  of  the  yield 
parameters  chosen  for  the  calculation  are  sunnarized  in  Table  IV.  The  values 
for  glass  were  obtained  from  shock  wave  compression  tests^3^  and  probably 
represent  compressive  fracture  rather  than  plastic  yield.  The  matrix  yield 
strength  was  obtained  from  the  values  measured  for  material  C (Table  II);  the 
bundle  yield  strengths  were  taken  to  be  the  same  as  the  matrix  in  the 
transverse  direction,  but  appreciably  larger  in  the  axial  direction  (large 
enough  that  kinking  always  preceded  macroscopic  plastic  yield,  in  accord  with 
observations) . 

Kinking  of  the  bundle  was  considered  to  Initiate  (Section  3.4) 
whenever  the  stress  differential  in  the  bundle,  <r22  - <rrr,  exceeded  a 
critical  value  of  400  MPa  (Fig.  5).  It  is  recognized  that,  for  statistical 
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reasons,  this  Is  likely  to  be  an  underestimate  of  the  kink  initiation  stress 
pertinent  to  the  Impact  of  a single  bundle  (i.e.,  an  overestimate  of  the  kink 
zone  size).  However,  It  is  a representative  value  suitable  for  present 
purposes.  Kink  propagation  was  assumed  to  occur  relatively  slowly  (because  it 
Involves  substantial  plastic  accommodation  of  the  fractured  fiber  segments  by 
the  matrix^);  hence,  kink  initiation  was  not  considered  to  produce  a 
significant  change  in  the  state  of  stress  within  a kinked  element  during  the 
penetration  phase  of  the  impact. 

A fracture  condition  was  also  invoked  to  ensure  that  the  tensile 
stress  did  not  attain  unacceptably  large  values.  A maximum  tensile  stress  cr 
was  attributed  to  each  element,  as  indicated  in  Table  IV.  The  element  was 
considered  to  fracture  when  the  average  tensile  pressure  exceeded  £.  This  is 
considered  satisfactory  for  the  present  problem  since  discrete  crack 
propagation  processes  are  not  involved. When  the  critical  tensile 
stress  in  an  element  has  been  exceeded,  the  stress  deviations  and  pressures 
are  set  to  zero.  Thereafter,  the  cell  can  still  support  compressive  stress. 
The  values  of  a selected  for  each  constituent  were  best  possible  values  from 
fracture  data  in  the  literature.  The  limitations  of  such  a selection  are 
recognized,  but  since  fracture  is  of  secondary  significance  in  the  present 
calculation,  the  choice  of  <r  is  not  critical  to  the  final  prediction  of  the 
kink  depth. 

A final  feature  of  the  calculation  was  the  introduction  of  an 
artifical  viscosity, which  was  required  to  maintain  numerical  stability  * 
in  regions  of  high  local  strain  rate.  This  introduces  an  error  into  the 
stress  calculation.  The  magnitude  of  the  error  is  most  significant  in  the 
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vicinity  of  the  shock  front.  The  error  appears  to  be  small,  as  discussed  in 
detail  elsewhere;^1*)  but  clearly,  the  influence  of  this  viscosity  term  must 
be  carefully  evaluated  in  each  case. 


5.2  RESULTS 

The  contact  stress  during  penetration  is  plotted  in  Fig.  19.  The 
value  approached  toward  the  final  penetration  stage  ( — 400  MPa)  is  reasonably 
close  to  the  measured  penetration  pressures  for  material  C (Table  III), 
providing  some  credence  to  the  material  properties  selected  for  the 
calculation.  Typical  values  of  the  axial  and  transverse  stresses  for  the 
center  cell  at  fixed  locations  along  the  bundle  are  plotted  in  Fig.  20.  The 
oscillation  in  the  stress  is  considered  to  derive  from  transverse  Inertia 
effects,  typically  encountered  for  Pochhammer-Chree  waves  propagating  in 


cylindrical  bars. 


The  development  of  the  kinks  during  penetration  is 


indicated  on  Fig.  21.  The  colls  in  which  the  kink  initiation  stress  has  been 
exceeded  are  Indicated  by  a horizontal  line.  Inspection  of  the  figure 
Indicates  that  an  appreciable  zone  develops  toward  the  final  penetration  stage 
In  which  only  the  outer  two  cells  of  the  bundle  experience  the  kink  formation 
stress.  These  cells  experience  the  requisite  stress  for  a very  short 
duration.  It  is  unlikely  that  the  duration  is  large  enough  in  most  Instances 
to  activate  fiber  fracture  (cf.  the  incubation  time  for  crack 
activation^11)).  A best  estimate  of  the  final  kink  depth  $ below  the  crater 
is  thus  4.4  mm  (Fig.  21).  This  Is  In  the  same  range  as  the  measured  values  at 
6,000  ms-1  (Fig.  16).  The  results  are  thus  quite  encouraging.  However, 
addUional  Insights  concerning  the  incubation  time  for  kink  initiation  are 
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clearly  needed  before  fully  quantitative  damage  predictions  can  be 
undertaken.  Also,  a means  for  incorporating  the  statistical  character  of  the 
kink  initiation  stress  must  be  identified. 


6.0  SUMMARY  AND  IMPLICATIONS 

Studies  of  carbon  composites  subjected  to  quasi-static  penetration 
and  to  projectile  Impact  have  demonstrated  strong  effects  of  microstructure  on 
the  extent  of  the  subsurface  damage  (manifested  as  kinks  within  both 
longitudinal  and  lateral  bundles).  Quasi-static  analyses  of  bundle 
deformation  and  the  resultant  fiber  curvature  have  indicated  the  important 
role  of  pre-existent  bundle  inhomogeneities  (initial  fiber  curvature,  matrix 
voids)  in  the  kink  formation  process,  axial  bundle  compression  also  emerges 
from  the  analysis  as  a prerequsite  for  kink  formation.  The  analysis  permits  a 
self-consistent  correlation  to  be  obtained  between  the  measured  kink  formation 
stresses,  the  penetration  pressures  and  the  quasi -static  extent  of  the  kink 
damage  In  longitudinal  bundles.  It  also  enables  the  morphology  of  the 
observed  kinks  and  the  kink  zones  to  be  rationalized. 

The  kinking  that  occurs  under  dynamic  conditions  can  only  be 
qualitatively  Interpreted  using  available  analyses  of  the  dynamic  target 
deformation.  Numerical  methods  are  needed  to  assist  In  the  development  of 
dynamic  damage  models.  An  Initial  attempt  at  the  numerical  computation  of 
kink  formation  In  the  longitudinal  bundles  has  thus  been  instituted  with 
encouraging  results.  Much  additional  work  Is  required  to  develop  a working 
model  of  Impact  damage. 
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The  observations,  measurements,  and  analyses  of  kink  formation 
provide  several  Important  Implications  for  the  design  of  damage  resistant 
composites.  If  the  depth  of  damage  in  the  longitudinal  bundles  Is  the  damage 
parameter  of  primary  concern  (i.e.,  the  damage  in  the  lateral  bundles  is  of 
secondary  concern),  as  might  be  the  case  for  erosion  resistance,  then  the 
preferred  bundle  geometry  Is  one  that  contains  the  maximum  feasible  proportion 
of  longitudinal  bundles  (with  a high  fiber  content):  the  bundles  being 
disposed  on  as  fine  a scale  as  possible  to  avert  deleterious  consequences  of 
the  local  anisotropy.  However,  If  the  extent  of  damage  In  the  lateral  bundles 
Is  also  of  Importance,  a more  uniform  array  of  bundles  may  be  preferred. 
Additional  work  Is  needed  to  quantify  this  requirement.  In  creating  the 
appropriate  bundle  geometry,  It  should  be  ensured  that  the  properties  of  the 
mlcrostructural  constituents  are  also  optimized.  The  matrix  has  a primary 
role  in  the  kink  Initiation  process.  Both  analysis  and  experiment  indicate 
that  a matrix  with  a maximum  possible  yield  strength  is  generally  preferred. 
Attempts  to  improve  the  Inelastic  properties  of  the  graphite  matrix,  or  to 
find  a substitute  material,  are  thus  to  be  encouraged.  However,  it  should  be 
recognized  that  there  may  be  an  upper  limit  to  the  damage  resistance  that  can 
be  Imparted  by  matrix  modification,  due  to  either  the  incidence  of  another 
damage  mechanism  (e.g.,  matrix  fracture)  or  to  an  Increase  in  the  penetration 
pressure  that  more  than  counteracts  the  Increase  In  kink  resistance.  Such 
eventualities  should  be  appreciated  If  discouraging  results  are  to  be 
avoided.  Another  primary  kink  formation  parameter  is  the  distribution  of 
pre-existent  bundle  Inhomogeneltles.  A reduction  in  the  amplitude  of  the 
large  extreme  of  the  bundle  inhomogeneities  should  lead  to  considerable 
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Increases  In  the  kink  resistance;  and  attempts  to  address  this  Issue  should 
constitute  a primary  direction  for  materials  development  studies. 
(Additionally,  the  fiber  fracture  strength  can  exert  an  Influence.  It's 
specific  role  Is  presently  rather  nebulous,  because  the  Importance  of  the 
fiber  strength  Is  contingent  on  the  magnitude  of  the  other  constituent 
properties— particularly  the  matrix  yield  strength  and  the  magnitude  of  the 
Inhomoyeneltles.  Nevertheless,  the  fiber  fracture  strength  should  be  Included 
as  a property  that  may  correlate  with  the  kink  resistance.) 

Finally,  the  consistent  trend  with  microstructure  cf  the  extent  of 
the  kink  zone,  the  yield  stress  In  constrained  compression  and  the  penetration 
pressure  is  emphasized.  This  suggests  that  these  relatively  simple 
quasi-static  tests  can  be  used  to  provide  a measure  of  a material's  impact 
damage  resistance,  suitable  for  materials  development,  studies. 
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APPENDIX  I 

COMPRESSIVE  DISTORTION  OF  FIBER  BUNDLES 


The  elastic  buckling  of  Initially  straight  laminates  has  been 
analyzed  by  Schuerch^®)  and  by  Rosen^7^  using  a strain  energy  method. 

For  the  mode  of  buckling  that  requires  extension  of  the  matrix  between  the 
fibers,  the  critical  buckling  stress  is  large,  and  given  by 


2V. 


3/2 


°-c- 


3(1-Vf) 


(Al) 


where  Is  the  volume  fraction  of  fibers.  However,  when  buckling  only 
subjects  the  matrix  to  axial  shear,  the  critical  buckling  stress  Is  reduced  to 


^m 

2Tl+v7(l-Vf) 


(A2) 


However,  even  this  stress  (~8  GPa  for  a typical  bundle,  Effl  » 5 GPa  and 
Vf  ~ 0.8)  Is  considerably  larger  than  the  measured  kink  Initiation  stress 
(0.2-0. 5 GPa).  The  onset  of  compressive  fiber  distortion  can,  therefore,  be 
reasonably  supposed  to  occur  In  regions  of  pre-existent  curvature  at  much 
reduced  stress  levels. 

The  deflection  of  initially  curved  fibers  in  a bundle  subjected  to 
compressive  deformation  Is  resisted  by  the  surrounding  matrix.  The  magnitude 
of  this  resistance  depends  on  the  presence  of  Interface  disbonds  and  matrix 
voids,  as  well  as  the  compressibility,  shear  modulus,  and  shear  strength  of 
the  matrix.  Consider  the  fairly  typical  case,  depicted  In  Fig.  22,  of  a 
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matrix  void  coincident  with  the  plane  of  fiber  curvature.  The  deflection  of 
the  adjacent  fiber  Is  then  dictated  by  the  applied  force  P and  the  resistance 

Imparted  by  the  material  on  one  side  of  the  fiber.  This  resistance  Is  ! 

j 

Initially  supplied  by  the  elastic  stiffness  of  the  matrix,  as  augmented  by  the  ] 

flexural  stiffness  of  the  nearest  neighbor  fibers.  The  subsequent,  Inelastic  j 

resistance  Is  related  to  the  compressibility  and  shear  strength  of  the  I 

matrix.  Assuming  that  inelastic  matrix  deformation  dominates  at  the  1 

I 

relatively  large  fiber  distortions  pertinent  to  kink  formation,  the  fiber  < 

behavior  can  be  analyzed  by  examining  the  solution  for  an  Initially  curved  ? 

1 

l 

elastic  beam  subjected  to  end  forces  P and  supported  along  Its  length  by  a j 

uniform  stress  q (Fig.  22).  Let  the  initial  curvature  y by  represented  | 

byU8)  j 

y0  » y sin  ir  (x/1)  (A3)  ] 

j 

where  the  dimensions  x,  1,  and  y are  shown  on  Fig.  22.  When  the  load  P is  j 

applied,  the  deflection  y Is  related  to  the  bending  moment  M by; 

! 

M - Py  - + gx^b  (A4)  j 


(A5) 


where  b Is  the  beam  width.  The  bending  moment  Is  also  related  to  the 
deflection  y'  due  to  the  deformation  by; 

M ■ - El  (d2y'/dx2) 

where,  y'  « y - yQ.  Combining  Eqs.  (A3),  (A4),  and  (A5)  gives 
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n dfjl  + 3|J.  . + py*  + pj  sin  ir(x/l)  • 0 (A6) 

■ 

Differentiating  Eq.  (A6)  twice  with  respect  to  x then  gives 


This  differential  equation  has  a general  solution  of  the  form 

y'  • A slnkx  + B coskx  + Cx  + D + F s 1 n (rrx/1 ) - qbx2/2P,  (A8) 

where  kc  ■ p/EI  and  A,  B,  C,  0,  F are  constants  of  Integration.  Applying 
the  boundary  conditions 

y'  ■ d2y'/dx2  » 0 


at  x * 0 and  x « 1,  we  obtain; 

-B  - D « qb/k2P 

A ■ (qb/k2P)[(cos{kl)-l)/s1n(kl)] 

C » qbl/2P 

F • yk2/  (»/l)2  - k2 


The  final  equation  for  the  deflection  Is  thus; 


y « -J- sin  Ux/n  - 

(x/1)  -k*  Pk‘ 

» - 


s<»(Kx)*cos(kx) 


m 


The  first  term  Is  due  to  the  Initial  curvature  of  the  beam  and  the  second  term 
Is  due  to  the  matrix  resistance.  In  applying  this  expression,  it  should  be 
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recall  ad  that  It  relates  only  to  the  deflection  that  occurs  when  the  matrix 
resistance  Is  Inelastic:  It  cannot  be  used  to  predict  the  Initial  beam 
deflections  (which  elicit  an  elastic  matrix  response). 

The  fiber  curvature,  *,  the  Inverse  of  the  radius  of  curvature  can 
now  be  derived  from  Eq.  (A9),  as, 


s1n(irx/l)  - 


(AlO) 


This  has  a maximum  negative  value  at  x * 1/2  given  by. 


A 

K « 


-y(V])4. 

. [(*/!)-*] 


+ 


qbf  1 - cos  (kl/2)  ' 

P cos  (kl/2) 


(All) 


The  stress  q that  restricts  the  development  of  curvature  In  initially 
distorted  fibers  transmits  lateral  forces  to  the  adjacent,  Initially  straight 
fibers  (Fig.  22).  This  Induces  premature  buckling  of  the  straight  fibers,  and 
propagates  the  zone  of  fiber  curvature.  The  displacement  of  the  force  P 
needed  to  develop  the  curvature  anticipated  by  Eq.  (AlO)  Is  not,  therefore, 
restricted  by  the  high  compressive  stiffness  of  the  surrounding  initially 
straight  fibers.  We  thus  propose  to  use  Eq.  (AlO)  as  a direct  measure  of  the 
relation  between  the  fiber  curvature,  the  applied  stress,  and  the  matrix 
resistance.  The  applied  stress  p Is  Introduced  by  assuming  that  the  stiffness 
of  the  fibers  still  greatly  exceeds  the  matrix  stiffness;  then,  for  a fiber  of 
square  cross  section^ 

P -pb2/Af  (A12) 
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where  Is  the  area  fraction  of  fibers  within  a bundle.  Noting  that  the 
moment  of  Inertia  Is  bA/l2,  the  maximum  curvature  becomes, 

A tt4£ (q/p) Af  1 - cos  [(1/b)  /3Vf  (p/Ef)  /Af 

* fl27r^l2Vf(p/Ef)(1/b)2/AfJ  5 cos  pl/b)  v/3Vf(p/Ef| 


( A13) 


Expanding  the  cosine  terms, 


cosx  « 1 - x2/2  + xA/24-... 


we  note  that  for  small  (l/b)z(p/Ef )/Af , 


-4a7  2 


3(q/Ef)( 1/b)' 


[>2-12(p/Ef)(1/b)2/Afj  2b  [l-3(l/b)2(p/Ef)/2Afj 


(A14) 


Rearranging,  enables  the  critical  applied  pressure  p for  fiber  fracture  or 

v» 

kinking  to  be  approximately  expressed  as, 


_ „ 2 /b\2  F A , 1 /VJv  1.512q 

EfAf 


(A15) 


where  *c  is  the  critical  (negative)  local  curvature  required  to  induce  fiber 
fracture  or  kinking.  For  small  matrix  shear  strengths  (q<<  Ef),  the 
Initial  dimensions  of  the  Inhomogeneity  (P,  | ) become  the  dominant  parameters, 


yielding, 


bl  ,2 


Pc'"  J Ef  Af  ^ 1 • ft 


(A16) 
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Indicating  that  p decreases  as  the  dimensions  of  the  inhomogeneity 

C 

Increase.  Relating  *c  to  a nominal  fiber  fracture  stress,  , , through; 


<r.  EfbK  rt 


(A17) 


the  critical  applied  stress  for  fiber  fracture  becomes, 


2 a r b 
7 AfEf  TT 


m 


(Ais) 


indicating  that  fracture  Is  Impeded  by  a large  fiber  fraction  and  fiber 
fracture  stress;  but  could  be  either  adversely  or  beneficially  influenced  by  a 
large  modulus  or  a large  fiber  radius,  depending  on  the  relative  magnitude  of 
the  other  variables. 

Alternatively,  relating  k to  a kink  formation  stress  (Appendix  II, 

V# 

Eq.  A22) , the  critical  uniaxial  compressive  stress  for  kinking  becomes, 


(A19) 


_k  2 « c b .2  w-yD  . 1 f/urn 

pc  * j AfEf  pr  ‘ - 5^-  Vvf  ( 1-v,  ■) 


In  this  instance,  a large  matrix  yield  strengths  replaces  the  fiber  fracture 
strength  as  a dominant  parameter:  while  the  other  parameters  (y,  1,  b,  E,) 
occupy  similar  roles. 
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APPENDIX  II 
KINK  INITIATION 


The  Initiation  of  a kink  Is  expected  to  ensue  whenever  a series  of 
correlated  fiber  fractures,  that  satisfy  the  condition  for  mechanical 
instability,  are  developed;^  usually,  this  will  not  be  coincident  with  the 
first  fiber  fracture. ^ Consequently,  the  mechanics  of  kink  formation  are 
intricate,  and  are  too  complex  to  attempt  to  derive  definitive  relations  at 
this  juncture.  Some  Important  insights  can  be  gained,  however,  by  examining 
the  fracture  of  laminates  consisting  of  a brittle  elastic  and  a ductile 
constituent.^9^  For  fiber  fracture  represented  by  the  relation 


*(S)  » 1 - exp  CS/SQ)m  dA  ] (A20) 

where  , m,  and  SQ  are  distribution  parameters  (cf.  Eq.  (2))  and  A Is  the 
surface  area  of  stressed  fiber,  the  critical  tensile  stress  c at 
instability  Is  given  by  the  relation/3,19^ 





(A21) 


where  w Is  the  width  of  the  bundle.  This  Is  a complex  relation  In  errand 
Insights  Into  the  factors  that  Influence  o£  can  only  be  obtained  for 
specific  situations.  An  Important  case  Is,  «Ef  and  m < 10.  For 
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example,  when  m ■ 5 (an  experimental  value  for  material  C is  8^), 
expanding  the  logarithmic  term  in  Eq.  (A21)  gives 


EfVf  (1  - 


m 


V.) 


( A22 ) 


In  this  case  the  fiber  fracture  parameters  ($Q,  A ) are  relatively 
unimportant:  the  matrix  yield  strength  is  the  dominant  Influence.  In  other 

i 

situations,  however,  the  fiber  fracture  parameters  can  exert  a more  important 
influence. 
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The  Intrafiber  stresses  that  derive  from  an  imposed  fiber  curvature 
are  dictated  primarily  by  the  elastic  properties  of  the  high  stiffness 
fibers.  The  role  of  the  matrix  is  to  maintain  a separation  between  fibers,  by 
virtue  of  their  resistance  to  shear  and  compression;  it  has  little  effect  on 
the  stresses  within  the  fibers  (i.e.,  a maximum  axial  shear  stress  given  by 
the  matrix  shear  strength,  r ).  However,  the  bundle  displacements  and  fiber 
curvatures  that  occur  in  response  to  an  imposed  force  can  be  strongly 
Influenced  by  the  properties  of  the  matrix.  These  properties  of  the  composite 
suggest  the  following  method  for  evaluating  the  Incidence  of  fiber  fracture 
and  kinking.  The  approximate  displacements  of  the  bundles  are  first 
determined  from  relations  derived  for  isotropic  elastic  or  elastic/plastic 
media  (most  pertinent  to  three-dimensional  composites);  allowing  the  fiber 
curvatures  that  develop  in  response  to  the  bundle  displacement  to  be 
obtained.  Then  the  tensile  stresses  within  the  fibers  are  deduced  from  the 
fiber  curvatures  and  equated  to  the  fiber  fracture  stress  or  the  kink 
Initiation  stress. 


1.0  BUNDLE  DISPLACEMENT 

The  medium  response  during  the  quasi-static  elastic/plastic 
penetration  of  a half  space  by  a rigid  sphere  has  been  shown  to  conform  quite 
closely  to  two  relatively  simple  limit  solutions. The  elastic/plastic 

34 

C/0726A/ebs 


Roskwi#  Internationa! 

Sconce  Center 


solution  for  a pressurized  spherical  cavity  provides  a reasonable  description 
close  to  the  projectile  Interface;  while,  the  response  in  the  elastic  zone 
remote  form  the  Interface  Is  most  effectively  described  by  the  elastic 
(Hertzian  and  Bousinesq)  solution  for  a half  space.  The  response  at 
intermediate  locations  can  be  estimated  by  interpolation.  The  radial 
displacement  within  the  elastic  zone  close  to  the  projectile,  predicted  by  the 
spherical  cavity  solution, is, 


rj  (1+0 


3<E  >r‘ 


(A23) 

(r  > rp) 


where  v Is  Poisson's  ratio,  < E > is  the  average  Young's  modulus  of  the 
composite,  r is  the  radius  of  the  plastic  zone,  <r  is  the  yield  strength 

r J 

and  r is  the  distance  from  the  center  of  the  cavity.  The  plastic  quantities 
r and  <r  are  related  to  the  conditions  that  exist  within  the  zone  of 

p y 

intense  fiber/matrix  distortion  near  the  projectile  interface.  Fig.  9. 

However,  we  note  that  rp  for  large  plastic  deformations  rapidly  attains  a 
limiting  value,  related  to  the  current  cavity  radius  r^,  by,^2^ 

(rp/rj)3  - <E>/3try  (1  - * ) (A24) 

Substituting  this  relation  for  rp  Into  Eq.  (A23),  the  radial  displacement 
becomes , 

(ur/r)  -i  (1  + *0/(1  - *0  (r^/r)3  (A25) 
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The  displacements  thus  depend  on  the  cavity  radius,  but  do  not  exhibit  a 
specific  dependence  on  the  plastic  properties  of  the  material.*  The  analogous 


large  deformation  result  for  the  plastic  zone  is: 


.(10) 


ur/r  - 1 - |[1  - (rj/r)3]  [1  - 2 (1  - 2*)(cy  < E >)  ] _1  } 


nl/3 


(A26) 


For  the  typical  condition,  cr^  < <<E>  , Eq.  (A26)  reduces  to: 


uf/r  - 1 


f1  - < V>31 


1/3 


(A27) 


Aqain,  therefore,  there  Is  no  specific  dependence  on  the  plastic  properties  of 
the  material.  The  displacements  close  to  the  projectile  interface  thus  assume 
a particularly  simple  form  (Fig.  23). 

The  elastic  displacements  remote  from  the  projectile  Interface— given 
by  the  half  space  solution— are  not  radially  symmetric,  but  have  a horizontal 


component  ux  and  a vertical  component  uz  given  by: 


.(22) 


u . U 2*111,+  , iQ 

x 2 < E > x 


poa 


2 2 “1/2  2 / 2 2\ *3/2  i 

z(xz  + Zz)  - 1 + xzz  (xz  + zz)  (l-2i>)  1 


poa 


uz  “ 7TT? 


(1 


+ v)  z2  (x2  + z2)  ^ + 2(1  - ,2)(x2  + z2) 


-1/2 


(A28) 


♦This  result  applies  strictly  to  Incompressible  plasticity.  The  equivalent 
result  for  compressible  deformation,  as  might  pertain  in  porous  materials, 
has  not  been  derived. 
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where  pQ  Is  the  contact  pressure  and  a Is  the  contact  radius.  The  remote 
displacements  thus  exhibit  direct  dependencies  on  the  plastic  properties, 
through  pQ,  and  on  the  projectile  radius,  through  a.  The  displacement 
characteristics  along  the  surface  and  along  the  penetration  axis,  predicted  by 
these  limit  solutions,  are  schematically  Indicated  In  Fig.  23.  Note  that 
there  Is  a region  of  zero  displacement  along  the  surface,  because  the  cavity 
and  half  space  solution  predict  displacements  ux  of  opposite  sign. 

2.0  FIBER  FRACTURE  AND  KINKING 

The  general  curvature  of  fibers  within  bundles  can  be  derived 
directly  from  the  above  displacement  relations.  Local  curvatures  (e.g., 
around  orthogonal  bundles  at  Inhomogeneities)  are  superimposed  on  the  general 
curvatures.  The  Important  trends  are  encompassed  by  examining  longitudinal 
bundles  along  the  penetration  axis,  lateral  bundles  just  beneath  the  surface 
and  longitudinal  bundles  beside  the  penetration  axis. 

2.1  AXIAL  LONGITUOINAL  BUNDLES 

The  longitudinal  bundles  along  the  penetration  axis  are  subjected  to 
axial  compression;  there  Is  no  general  curvature.  However,  a local  curvature 
associated  with  mlcrostructural  Inhomogeneities  (or  elastic  buckling)  can 
occur,  as  a direct  consequence  of  the  compressive  deformation  (Appendix  I). 

This  curvature  Is  related  to  the  magnitude  of  the  deformation,  the  matrix 
shear  strength,  the  fiber  modulus  and  volume  fraction,  and  the  dimensions  of 
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the  Inhomogeneity.  The  elastic  deformation  is  given  In  the  vicinity  of  the 
projectile  Interface  by: 


- - [2/9(1  - u)j  (rj/z)3 


(A29c) 


and  In  remote  regions  by: 

fz  • - 3pQa2/2<  E>z2  (A29b) 

The  equivalent  stresses  within  the  longitudinal  bundles  are, 

Pz  - - Eb  [2/9(1  - »)]  (rj/z)3  - /3(rj/z)3  (A30a) 

Pz  - - 3 (V<E>)  pQa2/2z2  . £*P0(a/z)2  (A30b) 


These  relations  should  permit  the  measured  kink  stresses  (Section  3.2)  and 
penetration  pressures  (Section  3.3)  to  be  correlated  with  the  depth  of  damage 
In  the  axial  longitudinal  bundles  (Section  3.3).  More  fundamental  predictions 
of  kinking  and  the  relation  to  microstructure  Involve  a combination  of 
Eq.  (A30)  with  fiber  fracture  relations  or  kink  formation  relations. 


2.2  LATERAL  BUNDLES 

The  lateral  bundles  are  subjected  to  a general  curvature  related  to 
the  displacement  u?.  Near  the  projectile  Interface,  the  curvature  K In  the 
elastic  zone  can  be  obtained  from  Eq.  (A23)  as: 


* 


3 -A 


z — rrrr- 

(xz  + z2)  (1 


(A31) 


38 

C/0726A/ebs 


Rockwell  International 

Sc  tone*  Center 


SC5076.5FR 


The  curvature  has  a minimum  (Fig.  24)  at  x » 0 given  by: 

*m1n  * - 3 (1  + ^ r^/z4  (1  - p)  (A32) 


and  maxima  at  x ■ + 3/4  2,  given  by: 

«max  - 0.85  (1  + p)  r^/z4  (1  - v)  (A33) 


The  curvature  In  the  plastic  zone,  derived  from  Eq.  (A27),  has  similar 
characteristics.  In  the  elastic  zone,  remote  from  the  projectile,  the 
curvature  derived  from  Eq.  (A28)  is. 


-2p0a‘  (1  + p) 


<E>  (x‘  + z £) 


777 


8x4  + 3x2z2  + p (2x4  - 9x2z2  - z4) 


(A34) 


which,  near  the  surface  (z  <<  x)  becomes: 


4pQa*  (1  + p) 

— : — 3 — 

< E > x 


(4  +p  ) 


(A35) 


Additional  local  curvature  of  the  lateral  bundles  can  occur  In  regions  where 
the  deflection  u2  exceeds  the  elastic  limit  of  the  matrix  Interstice.  The 
orthogonal  lateral  bundles  then  exhibit  a greater  resistance  to  distortion 
than  the  adjacent  Interstices,  and  perturb  the  deflection  of  the 
bundle^’3).  The  local  curvature  Is  derived  from  the  equation  for  the  local 
deflection,  as  a function  of  distance  x'  from  the  local  support^22) 

Nj,  -f(f)  flz(f)-9-4(f)zJ  (A36) 


j 
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and  Is  given  by, 


(A37) 


where  6 Is  the  relative  displacement  of  adjacent  out-of -plane  lateral  bundles 
and  i Is  the  spacing  between  bundles.  The  maximum  curvature,  at  x'  * 0 Is, 

“t'jZ  (A38) 


The  local  displacement  6 Is,  at  most,  equal  to  the  differential  In  the  general 
displacement  uz  at  two  locations  x and  x +1  along  the  bundle.  For  example, 
the  relative  displacement  of  near  surface  lateral  bundles  (small  z)  remote 
from  the  projectile  Is  (from  Eq.  A28) 

a^|pQa2  (l-i»2)/<E>x  (x+I)  (A39) 


giving. 


^12p0a2  {l-h/  <E > x (x+!)I 


(A40) 


The  lateral  bundles  are  also  subjected  to  tensile  and  compressive 
deformations.  In  the  elastic  zone  near  the  projectile  (Eq.  A25), 


e 


X 


ix2) 


m 


♦ *£) 


(A41) 


The  deformations  are  thus  positive  (tensile)  for  x/z  i 1/  and  negative 
(compressive)  for  x/z  ^ 1/  V2  (Fig.  24).  The  tensile  deformation  of  an 
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Initially  straight  fiber  does  not  produce  fiber  curvature;  whereas,  the 
tensile  deformation  of  an  initially  curved  fiber  tends  to  straighten  the  fiber 
and  eliminate  the  curvature.  Conversely,  compressive  deformation,  as  noted 
above,  tends  to  generate  curvature  at  local  Inhomogeneities. 

The  stresses  that  derive  from  the  general  flexural  curvature,  the 
local  curvature  and  the  axial  deformation  can  be  added  to  obtain  the  total 
stress,  c rxx.  In  the  tensile  zone,  o-xx  can  be  equated  to  the  fracture 
stress  to  estimate  the  Incidence  of  fracture.  In  the  compressive  zone,  It 
must  be  recognized  that  the  compressive  deformation  Is  augmented  at  local 
Inhomogeneities  (Appendix  I),  requiring  that  the  axial  term  be  coupled  with 
Eq.  (A22)  before  the  onset  of  kinking  can  be  predicted  by  equating  or^  to 
the  kink  Initiation  stress. 


2.3  NON-AXIAL  LONGITUDINAL  BUNDLES 

The  longitudinal  bundles  beside  the  crater  are  subject  to  curvature 
associated  with  the  displacement  ux.  However,  near  surface  relaxations 
limit  the  development  of  a general  curvature  in  the  Immediate  subsurface,  as 
well  as  restricting  the  zone  of  tensile  deformation.  Otherwise,  the  fiber 
curvatures  are  similar  to  those  anticipated  for  the  lateral  bundles. 
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TABLE  I 

PROPERTIES  OF  MATERIALS  SELECTED  FOR  PRESENT  STUDY 


Material 

Reinforcing 

Fibers 

Processing 

Sequence 

Unit  Cel! 
Dimensions  (mm) 

Relative 
Bundle  Areas 

X 

y 

z 

X 

z 

A 

Pitch  VS0022 

\ 

CVD/A240 

0.76 

0.76 

1.14 

1 



1 

1 

B 

Pitch  VS0022 

CV0/LS24 

0.76 

0.76 

1.14 

1 

l 

1 

C 

Thorne!  50 

CVD/PItch 

0.76 

0.76 

0.84 

2 

2 

3 

D 

Thorne!  50 

CVD/PItch 

0.76 

0.76 

0.76 

1 

1 

5 

E 

Thorne!  75 

Pitch/CVD/PItch 

1.02 

1.02 

0.56 

1 

1 

13 
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TABLE  II 

APPROXIMATE  ELASTIC  AND  PLASTIC 
PROPERTIES  OF  THE  COMPOSITE  CONSTITUENTS 


Material 

Test  Method 

Apparent 
Young's  Modulus 

(GPa) 

Apparent  Uniaxial 
Yield  Strength  of 
Matrix  (MPa) 

Bundle 

Matrix 

Axial 

Transverse 

A 

In-SItu 

220 

10 

3 

18 

Isolated 

190 

12 

4 

B 

In-SItu 

240 

12 

3 

20 

Isolated 

200 

10 

4 

C 

In-SItu 

200 

8 

5 

30 

Isolated 

180 

10 

7 

I 
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TABLE  IV 

MATERIAL  CONSTANTS  USED  TO  MODEL 
GLASS,  BUNDLE  AND  MATRIX 


Material 

Property 

Glass 

Bundle 

Matrix 

Equation  of 

State  Constants: 

p0 

2.5  g/cc 

1.87  g/cc 

1.87  g/cc 

A 

34  GPa 

25  GPa 

4.2  GPa 

B 

2 GPa 

0 

0 

Elastic 

Constants: 

vr$* 

0.22 

0.18 

0.18 

vrz* 

0.22 

0.18 

0.18 

"0r* 

0.22 

0.03 

0.18 

E ♦ 
rr 

62  GPa 

8.0  GPa 

8.0  GPa 

Grz 

25  GPa 

21  GPa 

3.4  GPa 

Yield 

Constants:** 

<r  0 

400  MPa 

(cry)zz  « 800  MPa 

17  MPa 

yi 

1 

°y 

Um 

12  GPa 

1.9  x 1010  ergs/g 

(<ry)rp  - 17  MPa 

0 

7 x 1010  ergs/g 

Tensile  Failure 
Constants: 

A 

®zz 

100  MPa 

2 GPa 

90  MPa 

A 

<rrr 

90  MPa 

A 

90  MPa 

°90 

♦Specified  only  for  anisotropic  material 

♦♦Yield  In  shear  for  Isotropic  materials.  Failed  materials  have  zero 
yield  strength. 
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A schematic  diagram  Indicating  kink  formation  In  (a)  longitudinal 
bundles  and  (b)  lateral  bundles. 

The  apparatus  used  for  constrained  compression  tests: 

(a)  schematic  of  load  application  to  sample,  (b)  an  overview  of  the 
test  fixture,  (c)  the  sample  and  compression  blocks. 

A kink  formed  just  above  the  yield  stress  in  a constrained 
compression  test. 

A typical  load/deflection  curve  obtained  In  a constrained 
compression  test. 

Kink  initiation  stresses  In  the  longitudinal  bundles  for  the  five 
materials  listed  In  Table  I. 

Kink  Initiation  probabilities  at  a constant  of  15  MPa  for  each  of 
the  five  materials. 

The  kink  strength  distribution  g(S|<)  plotted  as  a function  of  the 
strength  levek  S^. 

Damage  in  the  longitudinal  bundles  in  material  C created  by  a 3 mm 
diameter  penetration. 

Damage  In  the  lateral  bundles  In  material  C created  by  a 3 mm 
diameter  penetration. 

Scanning  electron  microscope  overview  of  the  crater  In  material  A 
at  velocities  of  4,000  and  6,000  ms-1. 

Kinking  In  longitudinal  bundles  In  material  C Impacted  at 
■6,000  ms-*:  (a)  section  through  Impact  center,  (b)  section 

adjacent  to  impact  center. 

Kinking  In  longitudinal  bundles  In  material  B on  sections  through 
Impact  center:  (a)  4,000  ms-1,  (b)  6,000  ms'l. 

Kinking  In  longitudinal  bundles  on  section  through  Impact  center  at 
6,000  ms"1:  (a)  material  D,  (b)  material  E. 

Kinking  In  longitudinal  bundles  In  materials  A and  B at  the 
extremities  of  the  kink  zone  Indicating:  (a)  an  unusual  mode  of 
kinking  at  the  corners  of  out-of -plane  lateral  bundles,  (b)  kink 
Initiation  at  a bundle/matrix  void,  (c)  kink  Initiation  at 
Intrabundle  voids. 
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Figure  15  Kinking  In  lateral  bundles  In  material  C impacted  at  6,000  ms-1. 

Figure  16  The  depth  of  kinking  in  the  longitudinal  bundles  for  1 mm  stress 

projectiles,  as  a function  of  the  impact  velocity. 

Figure  17  The  extent  x of  kinking  in  the  lateral  bundles. 

Figure  18  The  scheme  used  for  the  numerical  computation  of  kinking:  (a)  the 

cylindrical  bundle  Impacted  axially  by  the  glass  sphere,  (b)  the 
grid  used  for  the  computation. 

Figure  19  The  time  dependence  of  the  interface  stress  predicted  by  the 
calculation. 

Figure  20  The  time  dependence  of  the  stresses  developed  in  the  bundle  at  two 
distances  z below  the  surface. 

Figure  21  The  development  of  kinks  within  the  bundle  during  the  penetration 
process. 

Figure  22  A tylcal  bundle  Inhomogeneity  and  the  parameters  used  to  estimate 
the  fiber  curvature. 

Figure  23  The  displacement  in  the  longitudinal,  z,  and  lateral,  x,  directions 
predicted  by  the  elastic/plastic  pressurized  cavity  solution  for  an 
infinite  solid  (small  x,  z)  and  the  elastic  solution  for  a half 
space  (large  x,  z),  indicating  possible  interpolations. 

Figure  24  The  curvature  of  lateral  bundles  beneath  the  penetration  zone;  also 
shown  are  the  regions  of  axial  compression  and  tension. 
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KINK  INITIATION  STRESS,  Sk  (MPa) 
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